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Abstract. We review the inherent structure thermodynamical formalism and 
the formulation of an equation of state for liquids in equilibrium based on 
the (volume) derivatives of the statistical properties of the potential energy 
surface. We also show that, under the hypothesis that during aging the system 
explores states associated to equilibrium configurations, it is possible to generalize 
the proposed equation of state to out-of-equilibrium conditions. The proposed 
formulation is based on the introduction of one additional parameter which, in 
the chosen thermodynamic formalism, can be chosen as the local minima where 
the slowly relaxing out-of-equilibrium liquid is trapped. 



1. Introduction 

The possibility of a consistent description of the thermodynamics of equiUbrium and 
out-of-equiUbrium (glass) supercooled liquids has been and it is an important research 
line ||, ^, ^, ||, H, ^, H, |9[_ 1^ . In recent years, the inherent structure (IS) formahsm 



by Stillinger and Weber III | has significantly contributed to the understanding of the 
physics of supercooled liquids and appears to offer a powerful and simple approach 
for developing a thermodynamics of out-of-equilibrium (OOE) states. Indeed, on one 
side, the IS formalism provides a transparent way for writing the partition function 
in terms of the "base" of the local minima (inherent structures) of the underlying 
potential energy landscape (PEL). On the other side, state of the art computer 
simulations provide the possibility of a statistically complete sampling of the PEL 
explored in equilibrium conditions in a wide temperature range. The numerical 
analysis of configurations extracted from the canonical ensemble allows us to calculate 
the energy depth of the ISs explored during the dynamical evolution, to characterize 
the volume of their basins of attraction, and to give estimates of their degeneracy 
(configurational entropy) . Eventually it allows us to directly estimate the free energy 
of the system in terms of landscape properties. 

Recently j9| it has been shown that it is possible to write down a model equation 
of state (EOS) expressed only in terms of quantities describing the statistical 
properties of the PEL. The crucial step in this process is the evaluation/modellization 
of the volume dependence of the number of basins, of their energy distribution and 
volume. The landscape based PEL-EOS is able to predict the thermodynamics of the 
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system in equilibrium j^] and, even more interestingly, in out-of-equilibrium — 
when the aging system explores states related to equilibrium configurations. 

Although the calculations we have explicitely performed deal with the simulation 
of the Lewis and Wanstrom (LW) model for the fragile molecular liquid orthoterphenyl 
(OTP) our results are general. Here we review and discuss some general 

implications of our results on the understanding of the thermodynamics of supercooled 
liquids and glasses. 



2. The constant volume free energy 

In pioneering papers [0 , Stillinger and Weber have shown that the partition function 
Z(T) at constant volume V can be written as 

Z{T) = j dejs n{ejs) e'^^'^ e-^^^^^<^^'^'^\ (1) 

where /? = l/fc^T, ejs is the depth of the local potential energy minima (IS) of the 
PEL, fl{eis)deis is the number of potential energy minima with energy between e/s 
and e/5 + dejs, and Fmb{eis,T) describes the free energy of the system constrained 
in one of the basin of depth e/s, averaged over all basins of depth e/5. Following 
Stillinger and Weber , a basin is defined as the set of points in configuration space 
which lead to the same local minima under a steepest descent path. The power of 
this formulation relies in the fact that the procedures used to associate to each system 
configuration the corresponding IS are operationally well-defined through constant 
volume minimization techniques. Numerical evaluation of the density of states in the 
local minima allows us to calculate the harmonic contribution to the basin free energy. 
Starting from Eq. (|l|) , the free energy of the system can be written as 

F{T) = (e/s(T)) - T5co„/((e/s(T))) + F„6(T, (e/s(r))); (2) 

here (e/5(T)) is the solution to the saddle point approximation to Eq. (|l|), 
and Sconf = ln(ri(e/5)) is the configurational entropy. Fmb, the intrabasin 
vibrational free energy, is usually written in the harmonic approximation as Fyif, = 
kBT{'^f^^\u{f}'hu!i{eis)) , where uJi{eis) is the z-th normal mode frequency (z = 
1...M)) evaluated at the inherent structure, and h is the Planck constant. The sum of 
the logarithm of the normal modes frequencies describes the volume (via the curvature) 
of the basin of attraction of the IS in harmonic approximation. 

Computer simulation results and theoretical insight |l^ provide us valid 
models for the two crucial quantities 0(e/s) and Fmb, namely: 

n{eis)deis = e'^'' == deis (3) 

V Sttct^ 

Fy^bieis,T) = kBT[{a + heis) - kBTln{hf3)]. (4) 

The hypothesis of a Gaussian landscape is supported by the consideration that, 
if no correlation length diverges, e/s can be thought of as sum of the IS energy of 
several independent subsystems. In this case the central limit theorem suggests that, 
since the variance of the energy distribution in each of these independent subsystem 
is finite, a gaussian distribution will describe the distribution of e/s values We 
note that this hypothesis will break down at the very low energy tail, where differences 
between the gaussian distribution and the actual distribution become relevant. The 
second hypothesis ( J^fLi ln(/3?iwi(e/s)) — a + bejs ) is not crucial, but it is supported 
by the results of numerical studies. 
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Figure 1. T-dependence of (e/g) (a) and Sconf (b) in the approximation of 
Gaussian distribution of basin depths and ejg-independence of the basin volume. 
Curves for a = 1.0 and three different values of (t^/Bq are shown. 



Substituting in Eq. (g) and solving, one obtains |lq| : 

{eis{T)) = {Eo - ba^) - a^/ksT (5) 



-TS^onfiieisim = keT -aN^+ba'- ^ (6) 

F,^b{T,{eis{T))) =Fo{Eo,T)-kBTba\b + p). (7) 

Therefore, F(T, V) is expressed only in terms of proper combinations of the 
parameters a, Eg, cr, a, and b which are related to the statistical properties of the 
PEL and to a particular relation between volume and depth of the basins. We 
also note that from a plot of (e/5) vs 1/T it is possible to evaluate and Eq. A 
comparison between numerical data and Eq. ^ allows us to estimate a. In the case 
where all basins have the same volume (& = 0), Eqs. (||) and ^) simplify considerably, 
and in terms of scaled quantities one obtains 

&„./(fe(r)»/tB = a-5^£^. (9) 

Within the Gaussian approximation, the lowest e/5 value ck — characterized by 
(ck) = — is the Kauzman energy 



{eK(TK))lEo = l- ^2a^, (10) 
and it is reached at a Kauzman temperature Tk given by 



ksTk/Eo = Y — — — . (11) 

The behavior of (e/s) and Sconf {T) as a function of T, in reduced units, is shown in 
Fig. ||. We note on passing that recent works by Speedy Q and Sastry ||l^ have 
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Figure 2. Comparison between tlie different contributions to tlie pressure 
calculated according to tlie theory (solid lines) and by MD simulations (symbols) 
for the LW model: (a) Inherent structures contribution; (b) Vibrational 
contribution. The curves have been shifted by n X 20MPa to avoid overlaps; 
(c) Total pressure. Details on the calculation of these quantities can be found in 
Refs. I 0. 

attempted to correlate kinetic fragility to thermodynamic fragility suggesting 
that a and a are the statistical properties of the PEL which control the material 
fragility. 

3. Equilibrium equation of state 

The generalization of Eq. (||) to the volume dependent case requires the determination 
of the volume dependence of Eq. i.e., the formulation of an ansatz for the 
joint probability 'P{eis,V) to find a value of ejs at a given volume V. We follow 
the equivalent approach of fitting simultaneously the l.h.s. of Eqs. and 
determined by MD simulations at different volumes This procedure allows us 
to calculate directly the volume dependence of the parameters a, Eo, cr, a, and b 
introduced above. 

Substituting in Eq. (||) we obtain F{T, V) for the model considered, and the 
equation of state can be finally calculated via P{T, V) — —dvF{T, V) at T constant. 
From Eq. (^) it is immediately clear that P can be split into two contributions: a 
configurational part, Pconf, related to the change in the number and depth of available 
basins with V, and a vibrational part, Pvib, related to the change in the volume of the 
basin with V as 

P(T, V, eis) - PconfiV, ejs) + P^T, V, e/s)- (12) 

Fig. H shows the comparison among the MD estimates of the different contributions to 
the pressure (symbols) and the predictions of the above theory, for the case of the LW 
model. The excellent agreement between the two sets of data confirms the validity of 
the procedure introduced above which provides us an effective equation of state for the 
system under study based on the statistical properties of the landscape as expressed 
in Eqs. (|) and (|). 
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Figure 3. OOE simulation protocols: (a) Pressure evolution after a T-jump at 
constant volume per molecule V = 0.345 nm'^ from 480 K to 340 K; (b) Volume 
evolution after a P-jump at constant termjerature T = 320 K from 13.4 MPa to 
60.7 MPa. Details can be found in Ref. [Q. 



4. Out of equilibrium equation of state 

The possibility of a proper tliermodynamical description of out-of-equilibrium systems 
has been widely debated ||, ||, |l^, |o|. In particular it has been recognized that 
this should be possible by adding one or more history-dependent parameters to the 
equilibrium equation of state. The arguments discussed above allow us to go further 
in this direction jl^. Indeed the only hypothesis we made is the validity of Eqs. (||) 
and (Q). In all the out-of-equilibrium conditions where these two conditions are met, 
i.e., the system (gently) driven out of equilibrium explores states which are typical at 
equilibrium, the theory is expected to hold at the expenses of adding one parameter 
to the equilibrium EOS. Looking at Eqs. (^)(||) and the choice of the basin depth 
e/s as the additional parameter turns out to be very natural. 



To the extent of this extension the validity of Eq. (12) in OOE conditions is 
crucial, allowing us to link Pconf and P^b to e/s and V . If this is the case, the 
knowledge of e/s and V is sufficient to calculate both Pconf ^ Pvit and their sum P 
according to Eq. (|l2|). Similarly, the values of P, T and e/s are sufficient to predict 
V , since both Pconf {eis, V) and Pvib{eis,T, V) can be estimated as a function of V. 
The predicted V is the value for which Pconf {eis,V) + Pyib{eis^T,V) matches the 
external (fixed) pressure. 

In Fig. ^ we show the comparison among MD results and the predictions of the 
OOE equation of state for two different OOE protocols via computer simulation. In 
particular, we consider in Fig. |3|(a) the case of a T-jump at constant volume, and in 
Fig. ^(b) a P-jump at constant temperature. In the first case the dynamical evolution 
of e/s together with the (fixed) values of V and T allow us to predict the dynamical 
evolution of P; in the second one, the time dependence of e/s together with P and T 
allow us to predict the evolution of P. 

An interesting representation of the aging processes discussed above is the 
parametric plot in the Pis — eis plane. In Fig. || we show the path followed 
by the aging system for the two protocols discussed above (Figs. |(a) and (b)). 
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Figure 4. Paths of the aging process in the Pis — ^IS plane for the OOE protocols 
considered in Fig. ^ (a) T-jump at constant pressure; (b) P-jump at constant 
temperature. The arrows indicate the time evolution direction. Details can be 
found in Refs. jij . Panels (c) and (d) report the comparison between the "basin 
volume" during tne aging process with the "basin volume" of the corresponding 
basin (same ejs and same V) explored in equilibrium conditions. The basin 
volume is described by ^^i^j^ ln(a;i(e/s)). 



Panels (c) and (d) report the comparison between the basin volume described by 
the quantity X)f=i lii(/9?iwi(e/s)) during the aging process and the basin volume of 
the corresponding basin (same eis and V) explored in equilibrium conditions. In all 
cases the agreement between the calculated quantities and the theoretical prediction 
is quite good confirming the validity of our approach. 

Fig. H(b) is of particular interest, showing that in the OOE dynamics following 
a pressure jump can be recognized two different regimes. For times shorter than the 
barostat time constant (see Ref. |10 for details) the system responses to the external 
increase of pressure in a solid-like fashion, i.e., the PEL basins initially populated are 
only deformed by the volume change. Only for longer times, when the pressure has 
reached the equilibrium value, the system starts to age among basins different from 
the original ones. 



5. Conclusions 



In this paper we have reviewed some recent results on a general approach to the 
thermodynamics of equilibrium supercooled liquids and glasses |l^. We have 
discussed how it is possible to formulate an equilibrium equation of state in terms of 
quantities describing the statistical properties of the potential energy landscape. These 
findings allow us to better understand the nature of the different terms contributing 
to the total pressure of the system, and fill the gap usually found among experiments 
(usually performed at constant P) and computer simulations (usually performed at 
constant V). 
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The generality of the hypothesis we have introduced, allow us to generalize our 
approach to out of equilibrium conditions. In all the cases where the introduced 
hypotheses are met, i.e., the system ages among states typical at thermodynamical 
equilibrium, it is possible to write down an OOE equation of state at the expenses 
of the addition of one more parameter. This quantity can be naturally chosen as 
the depth of the explored inherent structures. The correctness of this generalization 
has been checked under several OOE conditions. Its limits of validity, under more 
severe OOE conditions where more then one additional parameter is needed for a 
consistent description (like in the so-called Kovacs memory experiments ^) 
is currently under investigation. We foresee the possibility that, under large variations 
of the temperature and/or pressure, different part of the system will age with different 
speed producing, as a net result, a material characterized by a distribution of ejs of 
the composing subsystems different from the equilibrium one, and/or a material for 
which the relation between basin volume and depth is different from equilibrium. 
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